Breakdown of superfluidity of a matter wave in a random environment 



M. Albert, 1 T. Paul, 2 N. Pavloff, 2 and P. Leboeuf 2 

1 Departement de Physique Theonque, Umversite de Geneve, CH-1211 Geneve, Switzerland 
Univ. Paris Sud, CNRS, Laboratoire de Physique Theorique et Modeles Statistiques, UMR8626, F-91405 Orsay 

We consider a guided Bose-Einstein matter wave flowing through a disordered potential. We 
determine the critical velocity at which superfluidity is broken and compute its statistical proper- 
ties. They are shown to be connected to extreme values of the random potential. Experimental 
implementations of this physics are discussed. 
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I. INTRODUCTION 

The simplest and most intuitive definition of super- 
fluidity (SF) is the ability to move without dissipation. 
According to a perturbative mechanism proposed by Lan- 
dau, superfluidity is broken when the velocity of the flow 
exceeds a critical value V C L at which it is energetically 
favorable to emit elementary excitations. Though this 
mechanism has been explicitly verified in 4 He [![, 3 He-B 
and Bose Einstein condensates (BECs) Q.many ex- 
periments in 4 He Q, 3 He-A @ and BECs [1, hav e 
shown that the actual critical velocity V c is generally 
lower than V C L due to the occurrence of phase slips. In 
this scenario, SF is protected by an energy barrier which 
may be overcome by fluctuations of thermal (as first sug- 
gested by Iordanskii in the context of liquid He-II § and 
by Little for superconductors M) or quantal origin (as 
more recently observed in 4 He [lOj, in superconducting 
nano- wires 11 1 and possibly also in BECs [Hj]), leading 
to what is called a resistive state in the physics of super- 
conductors. 

In what follows we address the problem of determining 
the critical velocity for the breakdown of SF of a mat- 
ter wave moving in a disordered potential. The matter 
wave beam is formed by a guided BEC in a quasi one- 
dimensional (ID) geometry. We assume zero tempera- 
ture; it is well known that in this case superconductiv- 
ity and superfluidity are not destroyed by weak disorder. 
This is Anderson's theorem for non-magnetic impurities 
in superconductors [13l ] ; similar results (with a different 
physical mechanism) hold for BECs [14]. In the latter 
case, the phase coherence of the system is preserved in 
ID in the presence of a weak disorder as demonstrated in 
the experiments reported in Refs. [ill HH- As a conse- 
quence, one may study a simple scenario for breakdown 
of superfluidity in disordered BECs where phase slips are 
neither thermally nor quantum mechanically nucleated 
but rather have a dynamical origin : the barrier dis- 
appears at a given critical velocity. This mechanism is 
standard in the absence of disorder (see, e.g., [TtJ and ref- 
erences therein) and the extraordinary control achieved 
in the domain of atomic vapor has even allowed a direct 
observation of the nonlinear excitations nucleated above 
the critical velocity V c [E Q- However, to our knowledge 
there is up to now only one clear experimental evidence 



of dynamical breakdown of SF and of finite critical veloc- 
ity in the presence of disorder, obtained by studying the 
damping of dipole oscillations in an elongated BEC [l8j]. 
In our fully ID well as in the dipole oscillation 

experiments [l8j . an important issue is to understand 
the out-of equilibrium solutions of a nonlinear continu- 
ous system in the presence of disorder. In this context, 
the phase diagram of the fluid flowing through a quasi 
ID disordered potential U(x) of finite extent L was re- 
cently studied in Refs. [19l - [2l| . Here, we concentrate on 
the SF part of this diagram and more specifically on the 
description of the breakdown of SF when the velocity (or 
the length L of the disordered region) increases. 

We study two different types of disordered potential 
with opposite characteristics. The first one is a smooth 
potential whose typical spatial scale of variation is large 
compared to the healing length of the condensate. In 
this case a local density approximation holds and a local 
Landau criterion can be applied [22]. This mechanism 
reconciles the Landau approach with the phase slip phe- 
nomenon, because it predicts that SF is broken when the 
local Landau velocity is reached by emission of nonlinear 
excitations (solitons in our ID case). The second type 
of disordered potential consists in a series of point-like 
impurities and thus has, contrarily to the previous type, 
strong fluctuations on small spatial scales. We show that 
in this case the criterion can be adapted yielding - as in 
the previous case - very good agreement with numerical 
simulations. In both cases we explicitly compute the sta- 
tistical properties of the critical velocity V c and show that 
they are closely related to the extreme value statistics of 
the disordered potential. We finally discuss experimental 
realizations of our models. 

The system considered is a weakly interacting BEC 
transversely confined by a harmonic potential of fre- 
quency u>±. For simplicity, the disordered potential U 
is supposed to depend on a single spatial variable - the 
coordinate x along the axial direction of the guide. A sta- 
tionary flow of the system is then accurately described 
by a ID order parameter ip(x) obeying the nonlinear 
Schrodinger equation [H, [24j 
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Here, n(x) = ^(x)! 2 is the condensate density per unit 
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of longitudinal length, u is the chemical potential and 
g = 2hijj^a v is the nonlinear parameter (a > is the 
3D s-wave scattering length). In the low density regime 
(an -C 1) the density profile in the transverse direction is 
Gaussian-shaped and v = 1, whereas v = 1/2 in the op- 
posite high density regime (an 3> I ) where the Thomas- 
Fermi approximation holds for the transverse degree of 
freedom 1251. 



II. SUPERFLUID FLOWS 

In addition to the density n(x), it is convenient to char- 
acterize the flow by its velocity v(x) — |j [arg (ip)] x - We 
assume that the disordered potential U(x) takes sizable 
values only over a region of finite length L. In this case 
a SF flow corresponds to a solution of Eq. (fTJ with con- 
stant density no and velocity V at ±oo. The chemical 
potential /x then reads: 



/i = ^mV 2 +gn v 



(2) 



In the following we denote the chemical potential of a 
BEC at rest as fia (u = gn^). 

In absence of external potential, the SF solution corre- 
sponds to n(x) = no and v(x) = V for all x, and one can 
show that it is stable under a weak perturbing potential 
provided V < cq where 
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is the sound velocity of the unperturbed condensate. The 
condition V < c is exactly the Landau criterion for SF 
because in a BEC the velocity is precisely the speed 
of sound. 



III. SLOWLY VARYING DISORDERED 
POTENTIALS 

In the case of a slowly varying potential, that is, when 
the typical length of spatial variations of U(x) is much 
larger than the healing length £ = hj ^/mfio of the fluid, 
one can devise a local density approximation for describ- 
ing stationary flows. In this scheme, the flow verifies 
current conservation and local equilibrium. This reads 
n(x)v(x) = C st — n V and 



fj, = -mv 2 (x) + gn u (x) + U(x) 



(4) 



From these considerations it is easy to see that the ve- 
locity v(x) is determined by a simple algebraic equation 
[H: U(x) = u G(v(x)) where G(v) = I - (V/vf + 
%(V 2 —v 2 )/cq. This equation admits a solution provided 
U(x)//j,q is lower than the maximum of G. If this con- 
dition is violated the flow no longer admits a stationary 
solution, SF is broken and the flow becomes dissipative 



(one can show that in this case the obstacle described by 
the potential U(x) experiences a finite drag [26]). The 
maximum of G is reached when V^Cq = [v (x)] 2+v , which 
precisely reads v(x) — c(x) where c(x) = \J vg n v (x) / m 
can be termed the local sound velocity [compare with 
(13])] : SF is broken when one reaches the local Landau 
criterion. 

In this approximation, it is clear that SF will first 
break down at the point x — x m where the potential 
is maximum: U(x m ) = max[[/(i)] = U m . The condition 
v(x m ) = c(x m ) yields an explicit relation between U m 
and V r : 
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Everything now boils down to a problem of extreme 
value: once the statistical properties of the maximum 
U m of U (x) over [0, L] are known, the distribution of 
the critical velocities V c can be obtained readily through 
Namely, if one denotes by & L (U m ) the probability 
distribution of U m and by Pl(Vc) the corresponding dis- 
tribution of the critical velocity V c one has 
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The distribution function of C/ m (and thus also that of 
V c ) depends on the size L of the disordered region for 
the simple reason that the longer the disordered region, 
the larger the probability of finding a large maximum of 
U. Hence it is clear that, on average, the critical velocity 
decreases with increasing sample length L. Note also 
that, in this picture, the critical velocity is related to the 
local fluctuations of the disorder and not to the details 
of the correlations. 

In order to calculate the distribution of C/ m , the first 
step consists in mapping the problem of finding the maxi- 
mum value of a correlated continuous function to a prob- 
lem of a set of N discrete uncorrelated variables. Ac- 
cording to |27| this mapping can be done if the corre- 
lation function of U(x), characterized by a correlation 
length £ C) decays faster than a logarithm and provided 
L/£ c 3> 1. From now on we assume that these require- 
ments are fulfilled, and study the extreme value statis- 
tics of a set of iV = jL/£ c uncorrelated random vari- 
ables {Ui,U2,—,Utf} distributed according to the prob- 
ability distribution of the disorder potential p(U). 7 is 
a parameter of order unity that depends on the correla- 
tion function and has to be determined numerically [28j . 
If we denote by f(U) and ^(f/ m ) the cumulative dis- 
tribution functions of U and U m , respectively, we have 
J£(Z7 m ) = [f(U m )] N . Taking the derivative of this ex- 
pression yields the probability distribution £^ L (U m ) and 
then the distribution of the critical velocity through ([5]) 
and ©. 

To estimate the frontier between the superfluid and 
the dissipative regimes, we calculate the median L(V C ) of 
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FIG. 1: (Color online) Transport of a quasi ID BEC with 
velocity V through a Gaussian correlated disordered potential 
of extension L (the parameters v, £ c and £ take here the values 
v = 1, £ c = 5£ and £ = O.l/^o). Dark region: time dependent 
flow; light gray (light blue online) region: SF stationary flow. 
The (yellow online) solid line displays the boundary between 
the two regions as predicted by Eq. (J7]). The inset displays 
the cumulative probability distribution of the critical velocity 
V c for samples of two different lengths. Staircase functions: 
numerical computations; dot-dashed lines: theory. 

this distribution and obtain 

L{Vc) = I' 1 In 1/2 
4 \nf(U m (V c )) ' 1 ' 

where U m (V c ) is given by (JSJ. 

In order to check the validity of our approach, we have 
numerically determined the critical velocity from time 
dependent simulations of the Gross-Pitaevskii equation. 
Starting from the ground state in the presence of dis- 
order at zero velocity, we have adiabatically accelerated 
the disordered potential until it reaches a velocity V. For 
each V and L we consider 80 realizations of the random 
potential and determine the fraction P s of stationary so- 
lutions. This quantity is plotted in Fig. Q] using a gray 
scale [dark, P s — ; light blue (gray), P s = 1] as a 
function of the normalized variables L/£ c and V/cq. The 
inset displays Fl(V c ), the cumulative probability distri- 
bution of the critical velocity. Very good agreement is 
observed in the expected limit of validity of our approach 
(L/£ c ^> 1). Figure Q] is drawn in the case of a Gaussian 
disorder, p(U) = exp(-[/ 2 /2£ 2 )/V27r£ 2 , with a Gaus- 
sian correlation function (numerically we find 7 ~ 0.8). 
We have also checked the accuracy of our predictions for 
other types of disorder of experimental interest (Lorentz- 
correlated disorder and speckle potential). 



p = N/L. A > is the dimensionless strength of a 
scatterer. In the presence of such a potential, the lo- 
cal Landau criterion is no longer applicable because the 
density is not smooth. However one can devise an ap- 
proach adapted to this particular case. One first remarks 
that each impurity repels the condensate, the density of 
which reaches its lowest local value at the position of the 
6 peak. The region in space where the decrease in den- 
sity is the largest will correspond to configurations where 
two (or more) scatterers lie very close to each other. SF 
will break down, therefore, at the point where the local 
concentration of 5 peaks is maximum. 

The smallest length scale for density modulations of 
the condensate is the healing length £. The conden- 
sate is therefore not sensitive to details of the disor- 
dered potential on scales smaller than £. One thus di- 
vides the disordered region into B — L/t; boxes, each of 
size £, and replaces the m, delta peaks present in box 
i by a single effective peak of strength m, A. The goal 
is then to determine the probability distribution of in- 
tensity of the strongest of the effective peaks, because 
it is at this peak that SF will first be broken. One 
thus needs to calculate the probability distribution of 
M = max {mi, 7712, tob}, where m i = AT. 

The probability of finding m peaks in an interval of 
length £ is tt(to) = Qp m (l - p) N ~ m , where p = £/L. 
In the limit of a wide disordered region (L — > 00), 
the product pN remaining constant, the binomial law 
can be approximated by a Poisson law of parameter 
( = pN: 7r(m) ~ e^^ C, m /m\ . In this limit the mi's 
are uncorrelated and the cumulative distribution of M 
is J[(M) = f(M) B where f(M) = T{M + 1,C)/M! is 
the cumulative distribution function associated with the 
Poisson law [r(x, £) is the incomplete gamma function]. 

One can relate the SF critical velocity to the strength 
A c = MA of the strongest effective peak by using the 
criterion for the critical velocity obtained in [23J for a 
single peak: A c = K(V c /c ), where 

K{z) = ^ {-8z 4 - 20z 2 + 1 + (1 + 8z 2 ) 3 / 2 } 172 . (8) 
The probability distribution of V c is then 

Pl(V c ) = K '^ c/co) & L (M = Ae(Vc)A) , (9) 
A c 

where & L (M) = d,^ L (M) /dM. Defining again the typical 
critical velocity as the median, the boundary of the SF 
region corresponds to 



IV. A SERIES OF <5 SCATTERERS 

Another commonly used model of disorder is a po- 
tential formed by a series of <5-like impurities: U(x) — 

-Vo£ J2i=i 3{ x — x i), where the Xj's are uncorrelated ran- 
dom variables distributed between and L with density 



In 1/2 



ln/(A c (V c )/A) 



(10) 



Figure [5] displays the stability phase diagram and the 
cumulative distribution of the critical velocity. Here also 
we find excellent agreement with the numerical results. 



4 




FIG. 2: (Color online) Same as Fig. Q]but for a random poten- 
tial formed by a sequence of uncorrelated delta peaks (A = 0.1, 
pi = 0.1). 

V. DISCUSSSION 

There are several experimental possibilities for study- 
ing the statistical properties of the boundary of the SF 
region. The Gaussian disorder with zero average may 
be implemented experimentally in the case of micro- 
fabricated circuits, where the atoms are magnetically 
guided over a chip [2j|. Roughness and disorder in the 
circuits induce fluctuations along the guide which are 
typically Lorentzian correlated, with a correlation length 
£ c which decreases with increasing distance between the 
guide and the chip [3(3]. However, the most common type 
of experimental disorder is the so called speckle poten- 
tial, generated by a laser beam passing through a dif- 
fusing plate [3l|. One of the most appropriate set ups 
seems to be the one used in Ref. Q, where the critical 
velocity of a trapped Bose-Einstein condensate has been 



probed by sweeping a laser beam through it. The crit- 
ical velocity was determined from measurements of the 
amount of excitations related to the emission of solitons 
and linear excitations. Similar studies could be done by 
sweeping the laser beam through a diffusive plate at con- 
stant velocity, thus creating a moving speckle potential. 
Finally, the statistical properties of the SF breaking may 
also be studied by simply adapting our calculations to 
the damping of dipole oscillations, along the lines of the 
recent experiments of Ref. [l8[ . 

In conclusion we have stressed the link between the SF 
critical velocity in presence of disorder and the statistics 
of extreme events. We have developed simple models in 
two opposite situations where the disorder is either very 
smooth or composed of point-like impurities. In both 
cases the agreement between numerical simulations and 
our analytical model is very good. We note that, be- 
cause of the mapping of the statistical properties of the 
critical velocities to that of extreme events of an uncor- 
related sequence, in the limit L — > oo the distribution of 
V c tends to one of the universal distributions of extreme 
value statistics 32]. A possible extension of this work 
could be the application to fermionic superfluids in the 
BCS regime. In that case, the Landau critical velocity is 
related to the pairing gap, and an equation similar to ([5]) 
can be explicitly written down [l7j . 
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